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In the covariant description of confinement, one expects the ghost correlations to be infrared 
enhanced. Assuming ghost dominance, the long-range behavior of gluon and ghost correlations 
in Landau gauge QCD is determined by one exponent k. The gluon propagator is infrared finite 
(vanishing) for n = 1/2 (k > 1/2) which is still under debate. Here, we study critical exponent 
and coupling for the infrared conformal behavior from the asymptotic form of the solutions to the 
Dyson- Schwinger equations in an ultraviolet finite expansion scheme. The value for k is directly 
related to the ghost-gluon vertex. Assuming that it is regular in the infrared, one obtains k ~ 0.595. 
This value maximizes the critical coupling a c (n), yielding a™ ax ~ (4tt/N c ) 0.709 ~ 2.97 for N c = 3. 
For larger n the vertex acquires an infrared singularity in the gluon momentum, smaller ones imply 
infrared singular ghost legs. Variations in a c remain within 5% from k = 0.5 to 0.7. Above this 
range, a c decreases more rapidly with a c — > + as k — > 1~ which sets the upper bound on k. 
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I. INTRODUCTION 

In gauge theories without Higgs mechanism, particles 
carrying the global charges of the gauge group cannot 
strictly be localized. Localized physical states are nec- 
essarily neutral in QED and colorless in QCD. The ex- 
tension to all gauge invariant and thus physical states 
is possible only with a mass gap in the physical world. 
Then, color-electric charge superselection sectors do not 
arise in QCD and one concludes confinement. 

The necessary conditions for this were formulated more 
than twenty years ago. In the next subsection we briefly 
recall these conditions, and how they constrain the in- 
frared behavior of ghost and gluon propagators in Lan- 
dau gauge QCD. Based on linear-covariant gauges, their 
derivation may not fully be divorced from perturbation 
theory. Their essence is quite generic and summarized 
in the Kugo-Ojima criterion which should apply in one 
way or another, whenever some form of BRS cohomol- 
ogy construction does for gauge theories. One way to- 
wards a non-perturbative definition of the Landau gauge 
is provided via stochastic quantization for which the 
full 5-dimensional BRS machinery is in the garage. The 
time-independent diffusion equation of this formulation is 
closely related to the Dyson-Schwinger equations (DSEs) 
in 4 dimensions as we describe next. Some of the nec- 
essary extensions, which have already been implemented 
in previous DSE studies of infrared exponents for other 
reasons, imply the Kugo-Ojima criterion. We summarize 
these studies, and how they are confirmed qualitatively 
in this way, at the end of the introduction. 

In Sec. II, we set up the DSE structures relevant for 
our present study. We summarize general properties 
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of the ghost-gluon vertex, most importantly its non- 
renormalization and ghost-antighost symmetry in Lan- 
dau gauge, which will be essential for infrared criti- 
cal exponents and coupling later. We then present the 
ultraviolet subtraction procedure with the special care 
necessary to make sure that it does not artificially af- 
fect the infrared. Some confusion arose recently con- 
cerning the relation between asymptotic infrared expan- 
sions and the renormalization group which we first clarify 
in Sec. III. We then review the non-perturbative def- 
inition of the running coupling that is based on the 
non-renormalization of the ghost-gluon vertex in Landau 
gauge, and show that in 4 dimensions it approaches a 
constant a c in the infrared whenever this vertex has an 
asymptotic conformal behavior also. As a byproduct of 
the vertex non-renormalization, the infrared behavior of 
both propagators thereby results to be determined by one 
unique exponent k in any given dimension. The general 
machinery to determine infrared critical exponent and 
coupling is outlined in Sec. IV. There, we also discuss the 
results with an additional regularity assumption on the 
vertex in the infrared, which in 4 dimensions leads to the 
values k « 0.595 and a c = a™ ax ~ (iir/N c ) 0.709 w 2.97 
for N c = 3. We furthermore discuss the infrared transver- 
sality of the vertex and show how this resolves an appar- 
ent contradiction with a previous study. 

We then discuss more general vertices involving an ad- 
ditional exponent which controls singularities in its ex- 
ternal momenta to discuss bounds on a c and k. Thereby 
we will find that values of k smaller than that for the 
regular vertex imply infrared divergences in ghost legs, 
whereas larger ones lead to an infrared divergence of the 
vertex in the gluon momentum. While the latter can only 
come together with an infrared vanishing gluon propaga- 
tor, which will always over-compensate this divergence, 
the former add to the infrared enhancement of ghost ex- 
changes. In particular, this would have to happen for 
an infrared finite gluon propagator (with k = 0.5) as 
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presently favored by lattice simulations. 

Our summary and conclusions are given in Sec. V, and 
we include 2 appendices which may provide the interested 
reader with some more technical details. 



A. The Kugo-Ojima confinement criterion 

Within the framework of BRS algebra, completeness 
of the nilpotent BRS-charge Qb, the generator of the 
BRS transformations, in a state space V of indefinite 
metric is assumed. The semi-definite physical subspace 
V phyB = Ker Qb is defined on the basis of this algebra by 
those states which are annihilated by the BRS charge Qb- 
Since Q% = 0, this subspace contains the space ImQ^ 
of so-called daughter states which are images of others, 
their parent states in V. A physical Hilbert space is then 
obtained as (the completion of) the covariant space of 
equivalence classes KerQs/ImQs, the BRS-cohomology 
of states in the kernel modulo those in the image of Qb, 
which is isomorphic to the space V s of BRS singlets. It is 
easy to see that the image is furthermore contained in the 
orthogonal complement of the kernel (given completeness 
they are identical). It follows that states in ImQs do not 
contribute to the inner product in V phyB . 

Completeness is thereby important in the proof of pos- 
itivity for physical states jl], ^ , because it assures the ab- 
sence of metric partners of BRS-singlets, so-called "sin- 
glet pairs" . With completeness, all states in V are either 
BRS singlets in V s or belong to so-called quartets which 
are metric-partner pairs of BRS-doublets (of parent with 
daughter states). And this then exhausts all possibili- 
ties. The generalization of the Gupta-Bleulcr condition 
on physical states, Qb\iP) = in V phys , eliminates half 
of these metric partners leaving unpaired states of zero 
norm which do not contribute to any observable. This 
essentially is the quartet mechanism: Just as in QED, 
one such quartet, the elementary quartet, is formed by 
the massless asymptotic states of longitudinal and time- 
like gluons together with ghosts and antighosts which are 
thus all unobservable. In contrast to QED, however, one 
expects the quartet mechanism also to apply to trans- 
verse gluon and quark states, as far as they exist asymp- 
totically. A violation of positivity for such states then 
entails that they have to be unobservable also. Increas- 
ing evidence for this has been seen in the transverse gluon 
correlations over the last years ||. 

But that is only one aspect of confinement in this de- 
scription. In particular, asymptotic transverse gluon and 
quark states may or may not exist in the indefinite metric 
space V. If either of them do, and the Kugo-Ojima crite- 
rion is realized (see below), they belong to unobservable 
quartets. Then, the BRS-transformations of their asymp- 
totic fields entail that they form these quartets together 
with ghost-gluon and/or ghost-quark bound states, re- 
spectively 0. It is furthermore crucial for confinement, 
however, to have a mass gap in transverse gluon correla- 
tions. The massless transverse gluon states of perturba- 



tion theory must not exist even though they would belong 
to quartets due to color antiscreening and superconver- 
gence in QCD for less than 10 quark flavors ||, ||, M. 

Confinement depends on the realization of the unfixed 
global gauge symmetries. The identification of gauge- 
invariant physical states, which are BRS-singlets, with 
color singlets is possible only if the charge of global gauge 
transformations is BRS-exact and unbroken. The suffi- 
cent conditions for this are provided by the Kugo-Ojima 
criterion: Considering the globally conserved current 

J«=d„F^ + {QB,D a J>c b } (with 0^ = 0), (1) 

one realizes that the first of its two terms corresponds 
to a coboundary with respect to the space-time exterior 
derivative while the second term is a BRS-coboundary. 
Denoting their charges by G a and N a , respectively, 

Q« = J dPx(diFSi + {Qb, D« b c b }) = G a + N a . (2) 

For the first term herein there are only two options, it is 
either ill-defined due to massless states in the spectrum 
of d v F^ vl or else it vanishes. 

In QED massless photon states contribute to the ana- 
logues of both currents in ([I]), and both charges on the 
r.h.s. in (^) are separately ill-defined. One can employ 
an arbitrariness in the definition of the generator of the 
global gauge transformations (||), however, to multiply 
the first term by a suitable constant so chosen that these 
massless contributions cancel. In this way one obtains 
a well-defined and unbroken global gauge charge which 
replaces the naive definition in (||) above ||. Roughly 
speaking, there are two independent structures in the 
globally conserved gauge currents in QED which both 
contain massless photon contributions. These can be 
combined to yield one well-defined charge as the gen- 
erator of global gauge transformations leaving any other 
combination spontaneously broken, such as the displace- 
ment symmetry which led to the identification of the pho- 
ton with the massless Goldstone boson of its spontaneous 
breaking § 0. 

If d v F® contains no massless discrete spectrum on the 
other hand, i.e., if there is no massless particle pole in 
the Fourier transform of transverse gluon correlations, 
then G a = 0. In particular, this is the case for channels 
containing massive vector fields in theories with Higgs 
mechanism, and it is expected to be also the case in any 
color channel for QCD with confinement for which it ac- 
tually represents one of the two conditions formulated by 
Kugo and Ojima. In both these situations one first has 
equally, however, 

Q a = N a = [Q B , J d 3 xDfc*} , (3) 

which is BRS-exact. The second of the two conditions 
for confinement is that this charge be well-defined in the 
whole of the indefinite metric space V. Together these 
conditions are sufficient to establish that all BRS-singlet 
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physical states are also color singlets, and that all col- 
ored states are thus subject to the quartet mechanism. 
The second condition thereby provides the essential dif- 
ference between Higgs mechanism and confinement. The 
operator D ab ~& determining the charge N a will in gen- 
eral contain a massless contribution from the elementary 
quartet due to the asymptotic field j a (x) in the antighost 
field, c a Xo —^°° jo. _| — , ^ m weak asymptotic limit), 

Here, the dynamical parameters u ab determine the contri- 
bution of the massless asymptotic state to the composite 
field g! ahc A^c h ^Z^ 00 u ab d^ b + - ■ ■ . These parameters 
can be obtained in the limit p 2 — > from the Euclidean 
correlation functions of this composite field, e.g., 

J d 4 x e*(*-») ( Dfc%x) gf bcd At(y)c c (y) ) = 

(^-^) U «V). (5) 

The theorem by Kugo and Ojima asserts that all Q a = 
N a are well-defined in the whole of V (and do not suffer 
from spontaneous breakdown), if and only if 

u ab ee u ab (0) = -S ab . (6) 

Then, the massless states from the elementary quartet 
do not contribute to the spectrum of the current in N a , 
and the equivalence between physical BRS-singlet states 
and color singlets is established pi 0, ^| . 

In contrast, if det(l + M) ^ 0, the global gauge symme- 
try generated by the charges Q a in eq. (|^) is spontaneu- 
osly broken in each channel in which the gauge potential 
contains an asymptotic massive vector field |2| . While 
this massive vector state results to be a BRS-singlet, the 
massless Goldstone boson states which usually occur in 
some components of the Higgs field, replace the third 
component of the vector field in the elementary quartet 
and are thus unphysical. Since the broken charges are 
BRS-exact, this hidden symmetry breaking is not directly 
observable in the physical Hilbert space. 

The different scenarios are classified according to the 
realization of the global gauge symmetry on the whole of 
the indefinite metric space of covariant gauge theories. If 
it is unbroken, as in QED and QCD, the first condition 
is crucial for confinement. Namely, it is then necessary 
to have a mass gap in the transverse gluon correlations, 
since otherwise one could in principle have non-local 
physical (BRS-singlet and thus gauge-invariant) states 
with color, just as one has gauge-invariant charged states 
in QED (e.g., the state of one electron alone in the world 
with its long-range Coulomb tail). Indeed, with unbroken 
global gauge invariance, QED and QCD have in com- 
mon that any gauge invariant localized state must be 
chargeless/colorless j2j. The question is the extension to 
non-local states as approximated by local ones. In QED 



this leads to the so-called charge superselection sectors 
||, and non-local physical states with charge arise. If in 
QCD, with unbroken global gauge symmetry and mass 
gap, every gauge-invariant state can be approximated by 
gauge- invariant localized ones (which are colorless), one 
concludes that every gauge-invariant (BRS-singlet) state 
must also be a color singlet. 



B. Infrared dominance of ghosts in Landau gauge 

The (2nd condition in the) Kugo-Ojima confinement 
criterion, u = — II leading to well-defined charges N a , can 
in Landau gauge be shown by standard arguments em- 
ploying Dyson-Schwinger equations (DSEs) and Slavnov- 
Taylor identities (STIs) to be equivalent to an infrared 
enhanced ghost propagator ||. In momentum space 
the non-perturbative ghost propagator of Landau gauge 
QCD is related to the form factor occurring in the corre- 
lations of Eq. (||) as follows, 

D oiP) = 4 * r a s . with U Q V) = S ab u(p 2 ) . (7) 

The Kugo-Ojima criterion, u(0) = —1, thus entails that 
the Landau gauge ghost propagator should be more sin- 
gular than a massless particle pole in the infrared. In- 
deed, there is quite compelling evidence for this exact 
infrared enhancement of ghosts in Landau gauge . For 
lattice calculations of the Landau gauge ghost propaga- 
tor, see Refs. [jio], |ll|, [l^]. The Kugo-Ojima confinement 
criterion was also tested on the lattice directly |0|. 

Lattice verifications of the positivity violations for 
transverse gluon states by now have a long history 
Q, [D| [L6|, |f7], [D§ Numerical extractions of their 
indefinite spectral density from lattice data are reported 
in p0[ |. As mentioned, however, this follows from color 
antiscreening and superconvergence in QCD already in 
perturbation theory |^| , and it is independent of con- 
finement. 

Its remaining dynamical aspect resides in the cluster 
decomposition property of local quantum field theory in 
this formulation [^|, |) . Within the indefinite inner prod- 
uct structure of covariant QCD it can be avoided for col- 
ored clusters, only without mass gap in the full indefinite 
space V. In fact, if the cluster decomposition property 
holds for a gauge-invariant product of two (almost lo- 
cal) fields, it can be shown that both fields are gauge- 
invariant (BRS-closed) themselves. With mass gap in 
the physical world, this then eliminates the possibility of 
scattering a physical asymptotic state into a color sin- 
glet consisting of widely separated colored clusters (the 
"behind-the-moon" problem) Q. 

The necessity for the absence of the massless particle 
pole in d v F^ v in the Kugo-Ojima criterion shows that 
the (unphysical) massless correlations to avoid the clus- 
ter decomposition property are not the transverse gluon 
correlations. An infrared suppressed propagator for the 
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transverse gluons in Landau gauge confirms this condi- 
tion. This holds equally well for the infrared vanishing 
propagator obtained from DSEs 21, |22|, ^3), and con- 
jectured in the studies of the implications of the Gribov 
horizon |25|, as for the infrared suppressed but pos- 
sibly finite ones extracted from improved lattice actions 
for quite large volumes [^6|, ^7], p8[ . 

An infrared finite gluon propagator with qualitative 
similarities in the transverse components appears to re- 
sult also in simulations using the Laplacian gauge [ p9[ . 
Related to the Landau gauge, this gauge fixing was pro- 
posed as an alternative for lattice studies in order to 
avoid Gribov copies |30| . For a perturbative formulation 
see Ref. |nj. Due to intrinsic non-localities, its renor- 
malizability could not be demonstrated so far. Devia- 
tions from the Landau gauge condition were observed 
already at 0(g 2 ) in the bare coupling in Ref. 1 32 . More- 
over, the gluon propagator was seen to develop a large 
longitudinal component in the non-perturbative regime 
p9j . In fact, compared to the transverse correlations, 
it seems to provide the dominant component in the in- 
frared. And it might in the end play a role analogous to 
that of the infrared enhanced ghost correlations in Lan- 
dau gauge. However, the precise relation with Landau 
gauge still seems somewhat unclear. It is certainly en- 
couraging nevertheless to first of all verify that no mass- 
less states contribute to the transverse gluon correlations 
of the Laplacian gauge either. 



C. Non-perturbative Landau gauge 



measure such that the expectation values of total deriva- 
tives w.r.t. the fields vanish, which is all we need to 
formally derive the same Dyson-Schwinger equations as 
those without restriction. 

In the stochastic formulation this restriction arises nat- 
urally because the probability distribution gets concen- 
trated on the (first) Gribov region as the Landau gauge 
is approached. Therefore, there should be no problem of 
principle with the existence of Gribov copies in the stan- 
dard DSEs. However, the distribution of the probability 
measure among the gauge orbits might be affected by 
neglecting (the non-conservative) part of the drift force. 
Ways to overcome this approximation are currently being 
investigated. Moreover, providing for a correct count- 
ing of gauge copies inside the Gribov region, the full 
stochastic equation will allow comparison with results 
from Monte-Carlo simulations using lattice implementa- 
tions of the Landau gauge in a much more direct and 
reliable way. In particular, this should be the case for 
the lattice analog of the stochastic gauge fixing used in 
simulations such as those of Refs. |l6|, |T§[ ]. 

Here, we restrict to the standard Landau gauge DSEs 
which are best justified non-perturbatively from the 
stochastic approach to be valid modulo the aforemen- 
tioned approximation. For their solutions, on the other 
hand, restricting the support of the Faddeev-Popov mea- 
sure to the interior of the Gribov region has the effect of 
additional boundary conditions to select certain solutions 
from the set of all possible ones which might contain oth- 
ers as well. Consider two invariant functions Z(k 2 ) and 
G(k 2 ) to parameterize the Landau gauge structure, 



A problem mentioned repeatedly already, which is left 
in the dark in the description of confinement within the 
covariant operator formulation presented so far, is the 
possible influence of Gribov copies |24| . 

Recently, renewed interest in stochastic quantisation 
arose, because it provides ways of gauge fixing in pres- 
ence of Gribov copies, at least in principle |3^, j34| . The 
relation to Dyson-Schwinger equations is provided by a 
time-independent version of the diffusion equation in this 
approach in which gauge-fixing is replaced by a glob- 
ally restoring drift-force tangent to gauge orbits in order 
to prevent the probability distribution from drifting off 
along gauge orbit directions. 

In particular, in the limit of the Landau gauge, it is the 
conservative part of this drift-force, the derivative w.r.t. 
transverse gluon-field components of the Faddeev-Popov 
action, which leads to the standard Dyson-Schwinger 
equations as clarified by Zwanziger |35|] . He furthermore 
points out that these equations are formally unchanged, 
if Gribov's original suggestion to restrict the Faddeev- 
Popov measure to what has become known as the in- 
terior of the first Gribov horizon is implemented. This 
is simply because the Faddeev-Popov measure vanishes 
there, and thus no boundary terms are introduced in the 
derivation of Dyson-Schwinger equations by this addi- 
tional restriction. Phrased otherwise, it still provides a 
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~k 2 



, D G (k) 



G{k 2 ) 
k 2 



,(8) 



in Euclidean momentum space of gluon and ghost prop- 
agator, respectively. When obtained as DSE solutions, 
suitable boundary conditions have to be satisfied by these 
functions, in addition. The following infrared bounds 
were derived by Zwanziger for each of the two as proper- 
ties of the propagators from the restricted measure: 

The observation that the "volume" of configuration 
space in the infinite-dimensional (thermodynamic) limit 
is contained in its surface lead to the so-called horizon 
condition which entails that the ghost propagator must 
be more singular than a massless particle pole in the in- 
frared [^, M EtI 



lim G _1 (fc 2 ) 

fc 2 ^0 



. 



(9) 



This condition is equivalent to the Kugo-Ojima criterion, 
u = — 1 for well-defined color charges in Landau gauge, 
c./., Eqs. (|) and (0) with G(k 2 ) = 1/(1 + u{k 2 )). 

From the proximity of the Gribov horizon in infrared 
directions Zwanziger furthermore concluded pq] that 



lim Z(k 2 )/k 2 = 



(10) 
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This removes the massless transverse gluon states of per- 
turbation theory as also required by the Kugo-Ojima cri- 
terion. The infrared vanishing of the gluon propagator is 
a stronger requirement than this, however. It currently 
remains an open question why this has not been seen 
in Monte-Carlo simulations as yet. An infrared suppres- 
sion of the gluon propagator itself, rather than Z(k 2 ), 
was observed for Landau gauge in f3§f| and, more con- 
siderably, at large volumes in SU (2) in the 3-dimenional 
case |39|, 40| as well as in Coulomb gauge Q. The 
3-dimensional results are interesting in that the large dis- 
tance gluon propagator measured there seems incompat- 
ible with a massive behavior at low momenta (that was 
noted also in |Il8|) . At very large volumes, it even be- 
comes negative |4C|, |ilf| . This is the same qualitative be- 
havior as obtained for the one-dimensional Fourier trans- 
form of the DSE results of [^2[ at small values for 
the remaining momentum c omp onents, c./., Fig. 4 of || 
vs. Figs. 2 of go) or 6 of @. Qualitatively, the dif- 
ferent dimensionality should not matter much here. On 
the other hand, the extrapolation of the zero momentum 
propagator in |4l| leads to a finite result which, however, 
still decreases (slowly) with the volume. This suggests 
the physical volumes may still be too small yet and that 
further study of the volume dependence of the zero mo- 
mentum gluon propagator might be necessary |71|. 



D. Infrared exponents in previous studies 

Within the standard BRS or Faddeev-Popov formula- 
tion, the functions in (|^) have been studied from Dyson- 
Schwinger equations (DSEs) for the propagators in var- 
ious truncation schemes of increasing levels of sophisti- 
cation JH Typically, the known structures in the 
3-point vertex functions, most importantly from their 
Slavnov- Taylor identities and exchange symmtries, have 
thereby been employed to establish closed systems of non- 
linear integral equations that are complete on the level 
of the gluon, ghost and quark propagators in Landau 
gauge. This is possible with systematically neglecting 
contributions from explicit 4-point vertices to the propa- 
gator DSEs as well as non-trivial 4-point scattering ker- 
nels in the constructions of the 3-point vertices || p2[ . 
Employing a one-dimensional approximation, numerical 
solutions were then obtained in Refs. [^2, 23 . 

Asymptotic expansions techniques were developed to 
analytically study the behavior of the solutions in the 
infrared. The leading infrared behavior was thereby de- 
termined by one unique exponent n « 0.92, 



z(k 2 ) k2 ~° 



with a renormalization group invariant a, see Sec. Ill A. 
The general bounds < K < 1 were established in Rcf. 
p2]| based on the additional requirement that Z and G 
have no zeros or poles along the positive real axis, i.e., in 



the Euclidean domain. Below, we will verify the pos- 
itivity for the leading infrared behavior of both these 
functions in the same range, independent of the one- 
dimensional approximation, and based on some few and 
quite generic properties of the ghost-gluon vertex alone. 

The infrared behavior in Eqs. ( pi] ) was later confirmed 
qualitatively by studies of further truncated DSEs. In 
Ref. [[l5| , a tree-level ghost gluon vertex was used in com- 
bination with a one-dimensional approximation which 
lead to a value of k w 0.77 for the infrared exponent of 
ghost and gluon propagation in Landau gauge. Then, in 
the first infrared asymptotic study of the ghost-gluon sys- 
tem without one-dimensional approximation, the value of 
k = 1 was obtained in Ref. ||4q| . There is, however, an 
issue about infrared transversality of the gluon propaga- 
tor, as we will explain below, which was not addressed in 
this study. As a result, the correct value for the tree- level 
vertex is the same as that derived herein for any ghost- 
gluon vertex with regular infrared limit, k 0.595, which 
was first reported for the tree-level vertex independently 
m Refs. H II- 

As we furthermore find in our present 
study, inconsistency arises for k — » 1 (from below), and 
this limit, the upper bound on k, is therefore excluded. 

With 1/2 < k, all these values of the infrared expo- 
nent share, however, the same qualitative infrared be- 
havior. The gluon propagator vanishes while the ghost 
propagator is infrared enhanced. Then, the Kugo-Ojima 
criterion and the boundary conditions (|)|), (|l0| ) are both 
satisfied. The horizon condition seems understandable, 
because the restriction to the Gribov region leads to a 
positive measure which was implicitly also assumed by 
requiring solutions without nodes in p2| . Here, depend- 
ing on the infrared behavior of the ghost-gluon vertex, we 
will find that this requirement could in principle be main- 
tained also for k < 1/2. Taken by itself, it only leads to 
< k < 1, and thus to the horizon condition (||). Even- 
tually, with decreasing k for values smaller than 1/2, in- 
frared singularities in ghost exchanges become too strong 
for a local field theory description. Around k = 1/2, how- 
ever, this argument is just not strong enough, and we can 
not turn it into an independent additional argument in 
favor of (|Io|) for an infrared vanishing gluon propagator. 



II. DYSON-SCHWINGER EQUATIONS 

The Dyson-Schwinger equations for the propagators 
of ghosts and gluons in the pure gauge theory without 
quarks are schematically represented by the diagrams 
shown in Fig. |l|. With infrared dominance of ghosts, the 
ghost loop represented by the diagram in the last line 
of Fig. [l] will provide the dominant contribution to the 
inverse gluon propagator on the l.h.s. in the infrared. 
In our infrared analysis we will concentrate on this con- 
tribution to the (renormalized) gluon DSE which reads 
in Euclidean momentum space with the notations of 0] 
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FIG. 1: Dyson-Schwinger equations for the ghost (top) and 
the gluon (bottom) propagator, diagrammatically. 



(color indices suppressed), 



D-i(fc) = Z 3 D\ v (k) - g 2 N c Z 1 x 

4 Gl(q,p)D G (p)D G ( q )G v (p,q) 



(12) 



(2^) 4 



where p = k + q, and the contributions from the four re- 
maining gluon loop-diagrams of Fig. |l| were not given ex- 
plicitly. D° is the tree-level propagator, D G is the ghost 
propagator, and G u is the fully dressed ghost-gluon ver- 
tex function with its tree-level counter part denoted by 
G ^. In the standard linear-covariant gauge the latter is 
given by the antighost-momentum, G®(q,p) = iq^. The 
DSE for the ghost propagator, without truncations at 
this point, formally reads, 



D G \k) = -Z 3 k 2 



+ g 2 N c Z 1 



d 4 q 

WT 4 



(13) 

ik^ D G (q) G v {q, k) D^ v (k - q) . 



The renormalized propagators, Dq and D^ v and the 
renormalized coupling g are defined from the respective 
bare quantities by introducing multiplicative renormal- 
ization constants, 



Z,D ( 



D 



G 



Z 3 D, 



Z g g = g ba 



(14) 



Furthermore, Z\ = Z g Z 3 Z 3 is the ghost-gluon vertex 
renormalization constant. Before we discuss the prop- 
erties of the ghost-gluon vertex, essentially the only un- 
known in ( |l2] ) and ([l]), we note the following: 

If we are allowed to assume that the leading contribu- 
tion to the inverse gluon propagator in (|l^) is completely 



determined by the ghost loop, this contribution must be 
transverse in Landau gauge. In other words, writing 

D~i{k) = eS^Zp^k 2 ) - k^k v Z R l {k 2 ) , (15) 

one should then have Zp(k 2 ) = Zp(k 2 ) = Z(k 2 ). Here, 
in particular, the leading infrared behavior as extracted 
from the ghost loop alone should not depend on whether 
we study Zp or Zr. With all other contributions sublead- 
ing, deviations from the transversality of the "vacuum 
polarization" , Zp = Zp in Landau gauge, should also be 
subleading. We will assess this by studying in parallel 
(we sometimes use D dimensions, normally D = 4 here), 



Z p \k 2 ) 



1 



D-l{k)V^{k) , 



Zp\k 2 ) 
with V^(k) = 6 ftp - 



(D-l)k 2 "* 



(16) 
(17) 



'''" kv , K„Ak; ,w n 1 ^ 



k 2 

respectively. Beyond the leading behavior as dominated 
by the infrared enhanced propagators within the ghost 
loop, there will in general be complicated cancellations 
between longitudinal contributions from various sources 
to ensure transversality of the gluons in Landau gauge. 
These sources can be due to the terms neglected, to trun- 
cations of vertices or to the regularization scheme em- 
ployed. The tadpole for example contributes only to Zp 1 , 
and so do the quadratic divergences with cutoff regular- 
ization. Beyond the leading order one therefore usually 
employed the 7?.-tensor in the contraction of the gluon 
DSE in most previous studies since that of Ref. pSfl . 

The tadpole contribution is a momentum-independent 
constant, so that it will necessarily be subleading as 
compared to the infrared singular ghost-loop, whenever 
that singularity is strong enough to lead to an infrared- 
vanishing gluon propagator, or Z(k 2 ) ~ (k 2 ) 2ft with 
K > 1/2 for k 2 -► 0. 

The infrared analysis we present below is independent 
of the regularization and both these reasons in favor of 
the 7?.-tensor do therefore not apply to our present study. 
Nevertheless, even with ghost dominance, exact transver- 
sality will in general only be obtained by including all dif- 
ferent structures possible in the ghost-gluon vertex that 
can contribute to the leading infrared behavior. 



A. The Ghost-Gluon Vertex in Landau gauge 

The ghost-gluon vertex is of particular importance in 
the analysis of the infrared behavior of the gluon and 
ghost propagators. We adopt the conventions of Ref. ||. 
The arguments of the ghost-gluon vertex denote in the 
following order the two outgoing momenta for gluon and 
ghost, and one incoming ghost momentum, c./., Fig. ^, 



Gf c (k,q,p) 

G a u c (q,p) 



(27r) 4 <T(fc + g 
gf abc G^q,p) . 



P )G a J> c (q,p), 



(18) 
(19) 
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M- > a 

§ k 



b y x c 

FIG. 2: Conventions for the ghost-gluon vertex G^ bc (k,q,p). 



Color structures other than the perturbative one assumed 
here were assessed for the pure Landau gauge theory on 
the lattice in Ref. Ell. In this study, there was no ev- 
idence seen for any significant contribution due to such 
structures which we will not consider henceforth. 

We parametrise the general structure of G M (q, p) which 
consists of two independent terms by the following form, 

G^q,p) = iq^A(k 2 ;p 2 ,q 2 ) + ik^B(k 2 ; p 2 , q 2 ) . (20) 

One might expect the second structure to be insignifi- 
cant in Landau gauge, since it is longitudinal in the gluon 
momentum k. This is not necessarily the case in Dyson- 
Schwinger equations, however, since the transversality of 
the vacuum polarization generally arises from cancella- 
tions of different longitudinal contributions as we dis- 
cussed above. 

For later reference, we recall two general properties of 
this vertex. The implications of these properties are ex- 
plored below. They both refer to the ghost-gluon vertex 
in Landau gauge: 

(N) Non-renormalization, Z x = Z g zl /2 Z 3 = 1 @, 
which entails that the vertex reduces to its tree- 
level form at all symmetric momentum points, 



G^(q,p) 



k- =q' 2 —p 2 



= G°Jq,p) , 



in a symmetric subtraction scheme. The gauge 
fields being purely transverse, however, there is a 
certain freedom left in the definition of the tree- 
level vertex. A priori, any form with r/ € [0, 1], 

Gl(Q,P) = VWn +ViPtt , and V + V = 1 . 

may be used equally for the Landau gauge. With- 
out further specification, for the functions in pi 
we first have, 



A(x;x,x) = 1, and B(x;x, x) 



V 



The condition on A is ^-independent. It expresses the 
essential aspect of non-renormalization and will be re- 
ferred to as (Nl). The condition on B depends on the 
ambiguity in defining the Landau gauge, as expected. We 
call this condition (N2). It reads B(x;x,x) — for the 



transverse limit of the linear-covariant gauge in standard 
Faddeev-Popov theory, as compared to B(x;x,x) = 1/2 
for the analogous limit of the ghost-antighost symmetric 
Curci- Ferrari gauge jM, see also Ref. H. These are the 
two special choices of particular interest, corresponding 
to r] = 1 and 77 = fj = 1/2, respectively. For renormaliz- 
ability and perturbative aspects of the latter, and for the 
geometry of the general 77-gauges, see Refs. fj52| . 

(S) The Ghost-antighost conjugation as part of the full 
Landau gauge symmetry, a semidirect product of 
SL(2, IR) and double BRS invariance, implies ||, 



A(x;y,z) = A(x;z,y) . 



(SI) 



While this holds for all 77, again, the B-function is 
more ambigous. It can not have definite symmetry 
properties in general. For the symmetrized Lan- 
dau gauge, based on the symmetric tree-level ver- 
tex with r\ = f\ = 1/2, the interactions with purely 
transverse gluons will preserve this exact symmetry 
of the Landau gauge, however. In the symmetric 
formulation we therefore expect to have an exactly 
ghost-antighost symmetric vertex also, 

G^(q,p) = G^ip.q) . 

Decomposing B = B + + B- with B±(x;z,y) = 
±B±(x; y, z), we then furthermore deduce, 



2B + (x;y,z) = A(x;y,z) 



(S2) 



In the fully ghost-antighost symmetric formulation we 
can thus express the vertex ( |20| ) in terms of the functions 
with definite (anti)symmetry as follows, 

G,(q,p) = ^±^A(k 2 ;p 2 ,q 2 )+ t k,B^k 2 ;p 2 ,q 2 ). 

(21) 

The B_ -structure is absent at tree- level and it vanishes 
at all symmetric points. Therefore, a symmetric vertex in 
converse to the logic above also requires r\ = fj = 1/2 for 
the tree-level vertex to be used as the in a symmetric 
subtraction scheme according to (N). 

At this point, it seems important to stress that, for 
the infrared exponent of Landau gauge QCD, only the 
infrared behavior of A is relevant. The critical exponents 
are of course independent of 77. As long as we concentrate 
on Z P via @ in the gluon DSE Q, the gluon legs of 
all vertices are transversely contracted. For all internal 
gluon lines this is automatically true by the transversal- 
ity of the propagator as in the ghost DSE ( |l3| ) for ex- 
ample, and for the external lines we just arranged it by 
hand. Thus, the 77-freedom in the tree-level vertex and 
the i?-structure of the full vertex are both irrelevant, as 
they should. In an infrared anaysis based on this mani- 
festly transverse system we might as well have standard 
Faddeev-Popov theory in mind with rj = 1, fj = 0. 

The only place where the rj-dependence and the B- 
structure do enter is the 7£-contracted gluon DSE. We 
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therefore introduce the generalized Landau gauge by the 
above modification of the tree-level vertex here as a 
purely technical tool to address the transversality issue, 
i.e., to compare Zp and Zp as obtained from Eqs. (|l6| ) 
and ([l?]), respectively. In fact, in order to reconcile ghost 
dominance with transversality, the result will be that for 
arbitrary values of 77 one must essentially have in the in- 
frared (indicated by the superscripts), 

5-(fc 2 ;pW) = ^A ir (k 2 ;p 2 ,q 2 ) , (22) 

BZ{k 2 ;p 2 ,q 2 ) = -t-^A*{k 2 ]P 2 ,q 2 ). (23) 

Inserting this into ( p0[ ) , small rearrangements reveal that 
the full ghost-gluon vertex therefore has to be transverse 
in the infrared itself, 



iq^pk - ip^gk 
k 2 



A* r (k 2 ;p 2 ,q 2 ). 



(24) 



This is in contrast to its perturbative limit where the 
.B-structure is suppressed. And it is now also indepen- 
dent of the choice of r\. Again, however, the transverse 
vertex is necessarily symmetric oc [iq^ + ip^) at a sym- 
metric point. In order to extend the subtraction scheme 
of (N) non-perturbatively into the infrared, and ensure 
transversality of the gluon propagator, we thus have to 
resort to the symmetric choice T] = 77 = 1/2 also (for 
which (|2|) follows trivially with (S2)). 



B. Truncated Slavnov- Taylor identity 

In pi| a Slavnov- Taylor identity of the standard linear- 
covariant gauge was derived to constrain the ghost-gluon 
vertex. Since the BRS transformations need some adjust- 
ments for other choices such as the ghost-antighost sym- 
metric gauges, as it stands this identity is valid for the 
case rj = 1 , fj = only. A generalization might be worth- 
while pursuing. However, considering the transversality 
of the full vertex in the infrared, this will not provide 
much additional information to be used in our present 
study, as we discuss in this subsection. 

Neglecting irreducible ghost-ghost scattering contribu- 
tions to the Slavnov- Taylor identity of Ref . |22| , and thus 
maintaining the disconnected contributions to the ghost 
4-point function only, a truncated Slavnov- Taylor iden- 
tity is obtained which, in terms of the two structures A, 
B in the vertex and the ghost propagator, reads, 



G(«)(^p 



A(y;x,z)-yB(y;x,z)) + (25) 

zG(x)G(y) 



G(z) 



Without the symmetry property (SI), a simple solution 
to Eq. (|2^) is given by 



A(x;y,z) 



G(x) 
G(z) 



This exact form was used for the ghost-gluon vertex in 
the study of Wilsonian flow equations for Yang- Mills the- 
ory in Ref. ]53]| . Implementing (SI) in addition, the most 
general solution to (]25|) can be written in the form, 

M \ G(X) G(X) rp 

A{x;y,z) = ——- + ——- - 1 + x j (x;y,z) , 
G(z) G{y) 

B(x;y,z) = ^-l + ^^f T (x;y,z). (27) 

The undetermined function f T thereby parametrises an 
unknown transverse contribution to the ghost-gluon ver- 
tex, kG T (q,p) = 0, of the typical type generally remain- 
ing unconstrained by the Slavnov- Taylor identities, 

G%q,p) = {iq„ kp - ip^ kq) f T (k 2 ;p 2 ,q 2 ) , (28) 

where k = p—q as before. We obtain, however, from ( [27] ) 
with (Nl) and (SI), respectively, 

f T {x; x, x) = , and f T (x;y,z) = f T (x;z,y) . (29) 

Since this function is otherwise arbitrary, in particular, 
in the infrared, the use of solutions (|27]) somehow seems 
less appealing for our present study which is concerned 
about the most general bounds on the infrared exponent 
that can be derived on the basis of as few and basic as- 
sumptions as currently possible. 

In the numerical solutions to the coupled system of 
truncated ghost-gluon DSEs presented in Refs. ]22|, p3[ , 
the form given in (^7|) was used for the ghost-gluon vertex 
with f T = 0. This solution to the truncated STI @ (for 
7] = 1) still satisfies (Nl), (N2) and (SI). Because it is 
not purely transverse in the infrared, it should be used in 
combination with the transversely projected DSE for Zp 
fro m (|l6| ), however, such that only the form of A(x; y, z) 
in (|27D matters. This causes ultraviolet problems in the 
numerical studies, see below. If the infrared transversal- 
ity of the vertex can be maintained on the other hand, 
by adding suitable transverse terms to a symmetric STI 
construction to satisfy (Q) above, for example, the 1Z- 
tensor may be used to contract the gluon DSE via jl7| ) 
by which these ultraviolet problems are avoided without 
doing harm to the infrared structure of the equations. 

We believe that this will be the way to proceed with 
the numerical studies of full solutions to truncated DSEs 
in future. In particular, this suggests to further develop 
the ghost-antighost symmetric formulation. 



C. Ultraviolet subtractions and infrared behavior 



With the parametrisation of the vertex in (20) we now 



B(x;y,z) = 0. (26) 



obtain for the ghost-loop contribution to the gluon DSE 
(O) the two alternative expressions from the contractions 
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according to Eqs. ( |l6| ) and (p7|), respectively, 



Z P {k 2 ) 



1 



Z R {k 2 ) 



Z:, 



= Z a 



, 9 2 iV c f A G(p 2 )G(q 2 ) 



k 2 p 2 q 2 
u2. „2 



5 2 iV c 



(2tt) 4 

</?>(% A(fc 2 ;<z 2 ,p 2 ) 

A G{p 2 )G(q 2 ) 



3 7 (2tt) 4 fcV<? 2 



(30) 



(31) 



[ V (qR(k)pA{k 2 ]q 2 ,p 2 ) - qK(k)kB(k 2 -q 2 ,p 2 )) + 
r)(^(fc)pA(fc 2 ;g 2 ,p 2 ) - pll(k)k B(k 2 ;q 2 ,p 2 ))} + ••• , 

where again p= k + q. One can see explicitly here that 
knowledge of both invariant functions is necessary for an 
infrared analysis based on the 7?.-tensor, c./., Eq. (31), 
while only the A-structure enters in Eq. (^) obtained 
with the transverse projector V . We furthermore allowed 
for the generalized tree-level vertex with r\ + f] = 1 dis- 
cussed in (N) which does not affect Eq. (|30|). This both 
makes the equation for Zp 1 particularly well suited for 
an infrared analysis, because then the invariant function 
A(k 2 ;p 2 , q 2 ), which parametrises the essential part of the 
vertex, is the only unknown remaining in the coupled sys- 
tem with the equally 77-independent ghost DSE (|l3|), 



G{k 2 



g 2 N c 



A Z{p 2 )G{q 2 ) 



(2tt) 4 k 2 p 2 q 2 

kV{p)kA{p 2 ;k 2 ,q 2 ) . 



(32) 



We used Z\ = \ for the Landau gauge in these equations. 
The ultraviolet divergences of the explicit loop-integrals 
are compensated by the renormalization constants Z3, Z% 
which we require to be infrared finite. For the ultraviolet 
subtractions, which are of course fc-independent, one then 
needs to make the following distinction: 

(i) In DSEs for propagators of massless or infrared 
enhanced degrees of freedom we can perform the limit 
k 2 — > 0. In the present case, we expect this for the ghosts, 
i.e., the l.h.s. of Eq. ( p2] ) will approach a finite constant 
Cg = linx^o G^ 1 (x) < 00 which can be zero, however. 
In such a case, the renormalization constant can easily 
be eliminated, e.g., here we then have, 



Z 3 = g 2 ^- D N t 
1 



D-l 
D 



d D q 
{2n) D 



(33) 



-Z( q 2 )G(q z )A(q 2 ;0,q 2 ) + C G 



We adopted the conventions of dimensional regulariza- 
tion with D — > 4~ here, but others may be employed 
as desired. If we assume regularity of A in the origin, 
we can conclude that A(x; 0,x) — > 1 for x — > from 
(Nl). Generally, if A(x; 0, x) does not vanish for x — > 0, 
then Eq. @ tells us that Z{q 2 )G{q 2 ) -> for q 2 -> 
in the infrared. Otherwise, Z% would be infrared diver- 
gent for D < 4. On the other hand, g 2 Z(q 2 )G 2 (q 2 ) is 
more and more recognized to be a good candidate for 



a non-perturbative defin ition of the running coupling in 
Landau gauge, see Sec. Ill A . However, this definition 
can only be reasonable, if we are able to arrange matters 
such that it does not vanish in the infrared (in fact, the 
running coupling must be monotonic to avoid a double 
valued /3-function with spurious zeros). 

So that with Z{q 2 )G 2 {q 2 ) ^ and Z(q 2 )G(q 2 ) -> we 
must have infrared enhancement of ghosts. In particular, 
Cg = in Eq. ( |33"|) and no such constant is then possible 
in a non-perturbative renormalization scheme. 

We can also reverse the logic here, and regard Cg = 
as an additional boundary condition on a set of possi- 
ble solutions to the DSE (|32"|). This then implements 
Zwanziger's horizon condition (||) to select the solution 
for the restricted Faddeev-Popov weight that vanishes 
outside the Gribov horizon |35fl , and that at the same 
time povides a positive definite G(k 2 ) > j22j. Once 
this selection is made, however, it then follows that 
Z(q 2 )G 2 (q 2 ) — * const, for q 2 — > (in D = 4 dimen- 
sions) which is a consequence of the non-rcnormalization 
of the vertex (Nl), as we will show in Sec. Ill C . So let us 
adopt C'g = and concentrate on the possible solutions 
with infrared enhanced ghosts from now on. 

With ([33]) in ([32j) we explicitly remove the ultravio- 
let divergence and obtain a manifestly finite equation to 
study the infrared behavior of G in D = 4 dimensions for 
a given form of A in the infrared, 



2 f A 3 Z(q 2 )G(g 2 ) 2 2 
9Nc Jj2^\4 A{q >°> q) 

~kV{p)k Z{ f 7 )G 7 { f ] A{p 2 -k 2 ,q 2 )) . [M) 



G(k 2 



k 2 p 2 q 2 

While the ultraviolet subtraction is rather simple with 
k 2 — > herein, without this subtraction, a naive infrared 
analysis will be aggravated by the ultraviolet divergences. 
Thus, the safe order of formal steps is to perform the 
ultraviolet subtraction before the infrared analysis in this 
case. The opposite order applies for the gluon DSE. 

(ii) In DSEs for propagators of massive degrees of free- 
dom or even infrared-vanishing correlations, the explicit 
ultraviolet subtraction is subleading in the infrared, and 
it cannot simply be extracted from the limit k 2 — * 0. 
This should certainly be the case for the transverse gluon 
correlations. The least we expect as a necessary condi- 
tion for confinement is the mass gap. The horizon con- 
dition implies an even stronger infrared singularity, as 
mentioned in the introduction. In either case we have 
Z _1 (fc 2 ) — > 00 for k 2 — > for the l.h.s. in DSEs such 
as ( |30| ) or (|3l]). The advantage is that the coefficients 
of the divergent terms in an asymptotic infrared expan- 
sion can be extracted without bothering about ultravio- 
let subtractions in the first place, since ultraviolet diver- 
gences will only occur at the subleading constant level 
in this expansion. To be specific, inside the ghost loop, 
the infrared enhanced asymptotic forms of ghost propa- 
gators (together the leading behavior of the ghost-gluon 
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vertex) will converge rapidly in the ultraviolet. Here the 
problem rather is to extract the necessary ultraviolet sub- 
traction without introducing by hand spurious infrared 
divergences. This problem had already to be dealt with 
in Refs. [^2], . What generally needs to be done in such 
a case is to reverse the order of (i) above, and to isolate 
the infrared divergent contributions to the gluon DSE on 
both sides prior to the ultraviolet subtraction, in order 
to define an infrared finite renormahzation constant Z%. 

Assume the coefficients in the infrared divergent terms 
of the asymptotic expansions for the ghost propagator 
and the vertex are known and the corresponding asymp- 
totic forms denoted by G lr (q 2 ) and A" lr (k 2 ;p 2 , q 2 ). We 
can then, e.g., in Eq. ( |30| ) subtract on both sides an ul- 
traviolet finite contribution of the form, 



1 



Z^(k 2 ) 



l q qV(k)q 



G' 



(2tt) 4 k 2 p 2 q 2 

"(p 2 )G lr (q 2 )A* r (k 2 ;q 2 ,p 2 ) 



(35) 



which we will actually use to determine Zp (k^) below, or 
an analogous expression for Z R r (k 2 ) from (pL). In both 
cases, this infrared subtraction in the gluon DSE has to 
be performed up to the necessary order such that 

! ' C A = const, for k 2 -> . (36) 



Z(k 2 ) Z ir {k 2 ) 

On the other hand, the fc-independent constant contribu- 
tion to the r.h.s. of the DSE is the one that contains the 
(overall) logarithmic divergence which is absorbed in the 
renormahzation constant Z3. To extract that from the in- 
frared subtracted gluon DSE, analogous to Eq. (^) from 
the ghost DSE in (i) above, we first rewrite the gluon 
DSE (pil), adding a zero via the definition in (|3F 

d D c, 



1 

= ^W) 

qV(k)q 



1 



D-l 



k 2 p 2 q 2 



7 ^{k 2 ) 
G(p 2 )G(q 2 )A(k 2 ;q 2 ,p 2 ) 



{2n) D 



(37) 



-G ir (p 2 )G ir (q 2 )A q 



V; 2 2 2 

(« ;q ,p . 



A further technical complication arises from the addi- 
tional power of q 2 /k 2 which prevents us from taking the 
limit k 2 — > explicitly herein. This problem is due to 
contributions which superficially contain quadratic ultra- 
violet divergences. Once the quadratic divergences have 
been taken care of, whenever this is necessary in a given 
scheme, we can then let k 2 — > for the logarithmically 
divergent contribution that defines the ultraviolet renor- 
mahzation constant Z3. For the 72-tensor, leading to Zr 
via Eq. (|l7| ) which is free from quadratically divergent 
contributions |48) , and for the generalized tree-level ver- 
tex G M = G° , as one special case to give an example 
without such complication, one readily verifies that 



/ 2(£>-2) 
\D(D+2) 

7 ( G ' V) 



rm 



{G ir (q 2 )Y 



(38) 



from k 2 — > in ( |3l| ) with the analogue of (35) for the 72- 
tensor. It is infrared finite by construction, and it gives 
the correct perturbative ghost-loop contribution to the 
gluon renormahzation constant Z3, ultraviolet divergent 
for D — > 4. Contrary to the ghost DSE renormahzation 
in Eq. |3^) , an additive constant is possible in (^) and is 
included in the terms not given explicitly here. It can be 
used to adjust the constant Ca in the gluon DSE which 
is subleading in the infrared, c./., (pf). 

To summarize, it should be possible to impose renor- 
mahzation conditions on the full propagators to equal 
the tree-level ones at an arbitrary (space-like) subtrac- 
tion point k 2 — p 2 > 0. To do this, we have two inde- 
pendent conditions at our disposal which fix the physi- 
cally insignificant overall factors in each of the two prop- 
agators. For massless tree-level propagators, we cannot 
without loss of generality extend this subtraction scheme 
to include k 2 = p 2 = 0, however, since that forced the 
full propagators to have the same massless single-particle 
singularity. Here, the necessity of a mass gap in the 
transverse gluon correlations entails that 1/Z(k 2 ) is in- 
frared divergent, and it is thus impossible to fix a multi- 
plicative factor by a subtraction at p = requiring that 
l/Z(0) be unity (or any other finite value). We can, how- 
ever, fix this factor by assigning the infrared subtracted 
l/Z(k 2 ) - l/Z lr {k 2 ) at k 2 = a certain value C A . 

That sets one of the two conditions available. Of 
course, the same argument applies to the ghost propa- 
gator at /j, = 0. In this case, it is because both, the 
Kugo-Ojima criterion and the horizon condition, tell us 
that the full ghost propagator should not have the sin- 
gularity structure of the free-massless tree-level one. In 
particular, with 1/G(k 2 ) — > we cannot fix the over- 
all factor with subtracting G" 1 at zero. As mentioned 
above, this case is different in that a non-vanishing con- 
stant contribution to the ghost DSE would be infrared 
dominant and cannot occur toghether with the infrared 
enhanced ghost correlations. To fix the multiplicative 
factor in the ghost propagator implicitly, we use the sec- 
ond of the two independent renormahzation conditions 
on the product of both propagators, 



GvmM 2 ) 



1 



(39) 



which can be used to define a non-perturbative running 
coupling in Landau gauge as we discuss next. 



III. RENORMALIZATION INDEPENDENT 
INFRARED ANALYSIS 

A. Infrared expansion and renormalization group 

Herein, we adopt the non-perturbative renormaliza- 
tion scheme introduced in Refs. |22|, ^3). To review this 
scheme briefly, recall that the formal solutions to the 
renormalization group equations for the gluon and the 
ghost propagator, e.g., for the latter this is Eq. (Al) in 
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App. |A|, can be written in the general forms, 



0(1) 



Z(k 2 ) = exp \ - 2 j" tk ' 9 dl ^ } f A (g(t k ,9)) , (40) 
G(k 2 ) = exp 4 - 2 £ (tfc ' 9 dZ ^ /G (^ (tfe; 5)) ^ (41) 



respectively. Here, = (lnfc 2 //z 2 )/2, and g(t,g) is the 
running coupling g(t,g), the solution of d/dtg(t, g) = 
(3(g) with g(0,g) = g and with the Callan-Symanzik 0- 
function, perturbatively /3(g) = -f3 g 3 + 0(g 5 ). The 
exponential factors are the multiplicative constants for 
finite renormalization group transformations (fx — > fjf), 

^(^)=exp{-2jf d*^}, (42) 

Z 3 ( M ',/i) = exp|-2£^^| , (43) 

with n' 2 = k 2 in (ftOl) , jil]), where 7,4(5) and 7g(<?) are 
the anomalous dimensions of the gluon and the ghost 
fields, respectively. 

The structure of (40), ( |4l| ) is summarized as follows: 
The momentum dependence of the propagator functions 
Z(k 2 ) and G(k 2 ) is completely determined by the running 
coupling evaluated at k 2 , which is renormalization group 
invariant, i.e., /(-independent, since (d/d\o.jj,)g(tk,g) — 
(fxd/dfi + (3(g)d/dg)g(t k ,g) = 0. We can therefore 
parametrise this momentum dependence by a function 
of the ratio of k 2 over a renormalization group invariant, 
dynamically generated momentum scale a cx A 2 CD . 

The ^(-dependence of the propagators, on the other 
hand, is then given only by the g = g(n) of the lower 
bound in the exponential renormalization factors. We 
can therefore always separate these two dependences, 
that on (g,fi) versus that on k 2 /a, in (40), pi] ) into 
multiplicative factors by conveniently choosing a go such 
that /1 2 = a at g — go, 



fj, 2 exp { 



dl 

W) 



(44) 



which, at the same time, defines a to be a renormalization 
group invariant momentum scale as promised. 

With this factorization of the (g, /^-dependence, we 
now make the Ansatz that the propagator functions Z 
and G have asymptotic infrared expansions in terms of 
k 2 1 a to some order TV, 

for k 2 1 a — > 0. Here we use a notation similar to that in- 
troduced in p4[. We note, however, that our expansion 



involves the RG invariant scale a while the renormal- 
ization scale fi was used in J54j. The difference can be 
absorbed in a redefinition of the coefficients e n , d n as we 
explain in App. Most importantly, this implies that 
our coefficients e„, d n are also (<?, /i)-independent. 

The non-renormalization (N) of the ghost-gluon vertex 

— 1/2 — 

in the infrared, c.f., Z\ = Z g Z 3 Z3 = 1 which was de- 
rived, in particular, for a symmetric subtraction scheme 
k 2 = p 2 = q 2 = /i 2 with fi 2 — > |]50| , now entails for the 
renormalization factors in Eqs. (|42|), (M3) that 



Zl /2 (ti',(i)Z3(fi',[i) = ! j 



* = (47) 



g 



in Landau gauge, which is equivalent to 
27c(ff) + 7a(s) : 



(48) 



This is, in fact, what allows to define a non- per turbative 
running coupling as introduced in Refs. E^, |23fl by, 



g 2 Z(^ 2 )G 2 ^' 2 ) = g' 2 =g 2 W/ii),g) 



(49) 



It reduces to the unique perturbative definition for large 
[i, fi' , is renormalization group invariant, dimensionless 
and thus as good as any non-perturbative definition can 
be. The fact that no constant of proportionality is in- 
volved in (|i^) implies a specific renormalization condi- 
tion. It corresponds to requiring the condition on the 
propagators, 

Z( M 2 ) = f A (g) , G(n 2 ) = f G (g) with f 2 f A = 1 , 

(50) 

which is incomplete, of course, to fix both their values 
separately at an arbitrary k 2 = /1 2 . The perturbative 
limits are, however, Ja,g ~~ * 1 > 5^0; corresponding to 
the perturbative MOM, 



Z(n 2 ) = 1 and G(n 2 ) = 1 



(51) 



for an asymptotically large subtraction point k 2 = [i 2 . 

With this so defined running coupling, by Eq. (fi"9|), the 
existence of an infrared fixed point, g(t, g) — > g c finite for 
t — >• — 00, then follows in Landau gauge to be in one-to- 
one correspondence with the scaling law for the leading 
infrared exponents of gluon and ghost propagation in the 
form (for D — A), 



e Q + 2S = 



(52) 



To make this explicit, consider the leading infrared be- 
havior from ([isl), ( |46l). with the exponential factors 
therein expressed by (|42j), (|43|) , for k 2 — > 0, 



Z(k 2 ) 
G(k 2 ) 



Zz(\[o,\i) e ( — 



Zz(\fB,\i) d — 
1 a 



So 



(53) 
(54) 
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which, from (|47|), entails that the infrared behavior of the 
running coupling a(k 2 ) = g 2 (tk, g)/(4ir) is given by, 



a(k 2 



?-Z(k 2 )G 2 (k 2 ) - 
4-7T 4l \ cr 



2 \ c o+2i5o 



(55) 



We therefore introduce a c = g^ody/^^) in the follow- 
ing. It represents the infrared fixed point, a(k 2 ) — > a c 
for k 2 — ► 0, which occurs exactly if (|5|) holds. 

The infrared scaling behavior in Eq. (52) was first ob- 
served in the solutions to truncated DSEs in Refs. |2^, 
p3f . It was recently derived from the ghost DSE in 
Ref . [glj . Therein, an additional assumption on the ver- 
tex was used which is a bit ad hoc and which is actually 
not necessary. In Sec. [II C below, we therefore present 



an alternative derivation of the infrared behavior (52) 
from Ref. {l7[ , which is based on the non-renormalization 
condition (Nl) for the vertex alone. 



B. Vertex Ansatz 

All we need in our infrared analysis is an Ansatz for 
the invariant function A which parametrises the relevant 
structure of the ghost-gluon vertex in Landau gauge. An- 
ticipating a conformal behavior in the infrared also for 
the vertex, we first write, 



A ir (k 2 ;p 2 ,q 2 ) 



(56) 



The non-renormalization condition (Nl) for the Landau- 
gauge vertex then leads to the constraint, 



I + m + n = 



(57) 



which will be implemented in our analyses, throughout. 
We will later also allow sums of terms of this form, in 
order to explore versions of this Ansatz which are sym- 
metrized w.r.t. the ghost legs. 

The scaling law for the infrared propagators in Eq. J52] ) 
then follows for such a sum of terms (|5^) via ([57]) only 
from (Nl), as we shall show in the next subsection. 

If we require, in addition, that A(k 2 ;p 2 , q 2 ) remain fi- 
nite when one of the ghost momenta vanishes, one of the 
terms of the form ( |56| ) must exist in the sum with either 



m = , I 



-n , or 



I = , m 



(58) 



All other possible terms must then vanish and thus have 
m > or I > 0, respectively. And if the finite contri- 
bution to A in that limit was to be in itself symmetric 
w.r.t. the two ghost momenta, one would only be left with 
A %T = 1 as in the tree- level vertex, since then I = m = 0. 
In a ghost-antighost symmetric sum of two terms on the 
other hand, we need one of each kind together with n < 
to avoid infrared divergent ghost legs. 

After the infrared scaling ([52]) and the general formu- 
lae for the infrared contributions of terms of the genuine 



form (p6|) , ( pTp will be derived, we will assume relations 
as in (pq), in addition. The joint infrared exponent K for 
ghosts and gluons then is a function of a single critical 
exponent n which is left as an open parameter in their 
vertex. To exemplify its influence, we will explicitly cal- 
culate this function for the following three Ansatze: 



(ii) 

(iii) 



A tr (k 2 ;p 2 : 



A^(k 2 ;p 2 , 



A* r (k 2 ;p 2 , 



k 2 

' n = 1 7 

1 

2 



fJ = n h? + -7 



(59) 
(60) 
(61) 



In the form (i) of ([39|) the Ansatz does not bother about 
ghost-antighost symmetry. For n = Sq = —k, this form 
contains the infrared behavior of the non-symmetric so- 
lution 
Ref. 
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(p6|) to the truncated STI ( (25| ) for the vertex of 
as a special case. In both symmetrized versions 
(ii) and (iii) one furthermore has n < 0, if infrared diver- 
gences associated with the ghost legs are to be avoided 
in the ghost-gluon vertex function. Version (ii) in ( |60| ) 
yields a finite and constant A(q 2 ; 0, q 2 ) = 1/2, while (iii) 
is an example for the possibility that A(q 2 ;0, q 2 ) = 0. 
Version (iii) in Eq. ( plj ) for n — <5o = — n includes the 
behavior of the symmetric solution (|27|) to the STI ( pq ) 
with f T = as obtained in Ref. p2[ ! All three versions 
satisfy (Nl), of course, and they all reduce to the tree- 
level vertex at n = 0. 



C. Unique infrared exponent from ghost DSE 

With the general form of our Ansatz ( p6| ) for the rele- 
vant part of the ghost-gluon vertex, we can now extract 
the leading contributions in the infrared on both sides of 
the ghost DSE @. Here, with Z ir and G lr denoting 
the leading infrared behavior of the propagators as given 
in (E33h and (p4j), respectively, we first note that 



2 j /2\ e o/2\ , 'o 



Z lr (p 2 )G lr (q 2 



. (62) 



where use has been made also of Eq. (pl7| ) for fi' = ^fo, 
g' = go- The l.h.s. of the ghost DSE (|34[) approaches, for 
k 2 -> 0, 



G- 1 (fc 2 )^Z 3 - 1 (^, / z)d - 1 



k 2 \- S " 



(63) 



To obtain the leading behavior at small k 2 /a of the r.h.s. 
in d34| ) , we replace the undetermined functions in the in- 
tegrand of (|34| ) by the from given in ( [62] ) and the Ansatz 
( jSq ) for the leading infrared behavior of the vertex. The 
difference in the integrand between the full functions, Z, 
G and A, and their asymptotic infrared forms, Z lr , G lr 
and A lr ', is subleading and it produces, upon integra- 
tion, terms that are also subleading in an expansion of 
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the r.h.s. of the DSE. This procedure is not restricted to 
the leading infrared behavior. It can straightforwardly 
be generalized for an infrared expansion to a given order, 
as long as all integrals in this expansion remain finite. 
This is true for the leading infrared forms as discussed 
above. When these are inserted, the integral in ( |34| ) con- 
verges in D = 4 dimensions. The role of the first term, 
however, just is to guarantee convergence of all necessary 
integrals at exactly D = 4. It is certainly necessary for 
the subtraction of the ultraviolet divergences when the 
full functions with their logarithmic momentum depen- 
dences are inserted. Here we only need this term, if we 
insist on a calculation in D = 4 involving convergent in- 
tegrals in every step. We can, however, obtain the same 
result by analytic continuation of integrals performed in 
D dimensions. Though this is of course precisely along 
the rules of dimensional regularization, this notion might 
be misleading here, since nothing is there to be regular- 
ized in the first place. It is adopted here for convenience 
only. Upon insertion of the leading infrared forms ( |6^ ) 
and (fitf), the first term in (|34|), in D dimensions, then 
yields a contribution to the ghost DSE proportional to 

which vanishes for any D in the analytic definition. 

Therefore, with (|62|) and in the integral for rele- 
vant part of the r.h.s. in D dimensions, and with (63) for 
the l.h.s., the leading contributions to both sides of the 



k' 2 



ghost DSE (| 

-So 



for k 2 



= -N r 4na n 



are readily extracted to yield, 

,2 \e +8o+Tr-2+l+rn+n 



(27T) 



D 



r 



kV( P )k 



S.-l+S +l 



X 



(65) 



1— £q— n 



Here, p = k ± q again, and we have 47ra c = g^egd 2 ] as 
introduced in ( |55| ) of Sec. Ill A. And just as for the usual 
replacement g — ► gfi 2 ~ D/ , the dimension of the coupling 
for D =f 4 has been taken care of in (j6^) by replacing 
9o ~ * <?o cr2 ~ D ^ 2 - The resulting explicit exponent of the 
scale can be combined with the corresponding exponent 
of an extra external momentum factor into {k 2 j o) D l 2 ~ 2 
which was added to the total exponent of this ratio on 
the r.h.s. of (|65[). The dimensionless integral therein is 
written in a form ready to apply the integration formula 
(Bl) of App. |B| The result is fc-independent as we will 
see explicitly in the next subsection. Therefore, from the 
non-renormalization of the vertex, condition (Nl) which 
implies I + m + n — 0, Eq. (^5|) for the leading infrared 
contributions to both sides of the ghost DSE entails that 



£o + 2<5 



L>/2 



(66) 



Thus, the infrared behavior of the running coupling from 
Eq. (pa) equivalently follows to be of the form, 



a(k 2 



2-D/2 



for k 2 







(67) 



which is infrared finite for D = 4. From now on, we there- 
fore parametrise the leading infrared exponents of ghost 
and gluon propagation by the joint exponent k again, in 
D dimensions, 



So 



-K, e = 2k + 2-D/2 



(68) 



Infrared enhancement of ghosts follows for exponents 
< k < (D — 2)/2 with the upper bound from the tem- 
peredness of local fields. In addition, (D — 2)/4 < k for 
the mass gap in transverse gluon correlations. Together, 
this then establishes the Kugo-Ojima confinement cri- 
terion in Landau gauge. Zwanziger's horizon condition 
furthermore excludes equality at the lower bound. We 
thus expect solutions for 



(D-2)/4<k< (D-2)/2 



(69) 



Once this exponent is determined, depending on the in- 
frared exponents of the ghost-gluon vertex, from Eq. (|65| ) 
we furthermore obtain, 



N c I ( G D) ( K ,l,n) 



(70) 



where Iq* (k,/, n) is a ratio of T-functions proportional 
to the integral in Eq. (Jog) which we determine next. 



IV. INFRARED EXPONENT FOR 
GHOST-GLUON SYSTEM, RESULTS 

The same procedure that led to (|6^) in the ghost case, 
can be applied to the gluon DSE. In this case, we simply 
insert the leading infrared behavior of the propagators 
from (||), (H) and the vertex @ directly into Eq. @. 
In D dimensions, this then analogously leads to 

(71) 



— J = N c 4na c ( — 



2 D TT D / 2 



where we used l + m + n = 0, and again, we conclude the 
relation for the exponents in Eq. (|66|). A comparison of 
Eqs. ( |70| ) and (|7l]) furthermore tells us that 



/^(K,/,n) = l¥}(K,l,n) 



(72) 



which determines the values allowed for the exponents 
n, I and k. The two dimensionless integrals 1^ and 
1^ for the transverse ghost-loop contribution are given 
byf 



I^ D) (K,;,n) 
2 D ir D / 2 



Iz D p ) ( K > l > n ) 

2 D TT D / 2 



(27T) 



kV(p)k 



-l — K+l 



(2nY 



qP{k)q 
{D-l)k 2 

l-K-n-l /, 2 
P 2 



-2k— n 

, (73) 



l+K-l 



(74) 
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To compute these, we first note that 



kV(p)k 1 / q 2 q 2 \ 1 (k 2 

p2 £2 



k 2 

qP{k)q 
k 2 



4 \p 



(£_ 

k 2 



El 

k 2 



2 2 
PI 



El + q 

k 2 k 2 p 2 

■A „4 



If q* p 
4 { + ¥ + fc 1 



With these in Eqs. ([73|) and (|74j), respectively, it is now 
straightforward to apply the formula in Eq. ( JB l| ) repeat- 



edly with suitable substitutions for the exponents a and 
(3. For each of the two integrals this leads to a sum of 
six ratios of F-functions with arguments that differ by 
integer values corresponding to the above six terms with 
different powers of the momenta in each of the two in- 
tegrals. After some tedious applications of the T's func- 
tional identity, these six ratios can be combined into a 
single one, to the effect that, 



J 



G 



(k, I, n) 



D - 1 



r( 



k + l)T(l + 2k + n)T{-K -l-n) 



r(f -2k- n)r(l +k- l)T{f + 1 + k + I + n) 



1 T(f -k + l)T(l - § + 2k + n)r(f - « - I - n) 

2 r(D - 2k - n)r(l + k - z)r(i + k + / + n) 

I 



(75) 
(76) 



Though these quite general results might look compli- 
cated at first, in fact, they are surprisingly simple. To 
appreciate this, consider the following special case first. 



A. Tree-level ghost-gluon vertex 

In this subsection we concentrate on the results in the 
event that the ghost-gluon vertex reduces to its tree-level 
form in the infrared, G% — G°. This corresponds to 
A tr = 1 and, for Zr, the generalized form with B" = fj. 
We thus set I = n = in the general result of Eqs. (75), 
(f76|), and omit the arguments I and n here. In D = 4 
dimensions, this amounts to 



4 4) («) 



4>) 



3 r 2 (-K)r(2«- 1) 
~2T(-2k)T(k- i)r(« + 3) 
1 T 2 (-k)T(2k- 1) 
2T 2 {k- l)r(4-2/c) ' 



(77) 
(78) 



The result in Eq. (77) then agrees with both versions 
given for the tree- level vertex in Ref. |4q |. There, another 
method was employed leading to a rather complicated 
form which consists of sums of confluent hypergeometric 
functions and which is quite difficult to simplify any fur- 
ther. We explain this method and the connection to ours 
in App. [B|, c./., Eq. (Bl) vs. Eq. (B7), in particular. The 
two different forms for the result in Ref. [ff6[ thereby arose 
from either choosing the internal gluon or the ghost mo- 
mentum as the integration variable in the loop. Here, 
both gives the same in the first place. The symme- 
try p 2 <-> q 2 in the integrand on the l.h.s. of Eq. (Bl) 
is manifest on the r.h.s by its explicit invariance under 
a -> a' = D/2 - f3 together with j3 -> /?' = D/2 - a. 

The selfconsistency condition (|72|), for both DSEs to 
yield the same value of the constant a c , is also imple- 
mented quite easily in Eqs. ( |77| ) and (|78|). In this case, 



for the tree-level vertex, one derives the condition, 

12(3-2k)(2k- 1) = (k + 2)(k+1) . (79) 



This is quadratic in k, and the two possibilities are, 



k = 



— (93 T Vl20l) « { 0.59535 . 1.3026}, 
98 V / 



(80) 



with one unique root in < K\ < 1 which we have un- 
derlined in (|80|). This result was first obtained indepen- 
dently in Rcfs. |35], ^7). The corresponding value of a c is 
given by, 



47T 



with4 4) ( Kl ) =I£(ki) 
last factor than the value of a 



2.9717 , for N c = 3 



(81) 



1.4096. It is smaller by this 
4ir/N c derived from 



Iq\k) = 1 for K = 1 with the tree-level vertex (and 
also D = 4) in Ref. [|46) . This disagreement is due to the 
difference between the transverse ghost-loop contribution 
Zp employed here, and the ^-contracted Zr in Ref. jl6). 
Our calculation for Zr, see below, with A = 1 and B = fj 
for the gerneralized tree-level vertex leads to 



r (4) 



%) = (4.-2) (l - 2^^-^) 4>) . (82) 

For rj = 1 or this is equivalent to the sum of twelve 
confluent hypergeometric functions given as the result 
for the same integral in Ref. [ fhS| . The D = 4 results 
for the tree-level ghost-gluon vertex are summarized in 
Fig. ||. In standard Faddeev-Popov theory r\ = 1, and 

transversality of the ghost loop, 4= = 4p> required 
k = 3/4. In order to tune Z R for transversality at the 
selfconsistent value of the tree- level exponent K\, on the 
other hand, we would need r\ sa 1.16 or —0.16 which 
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uon vertex, c./., Eqs. ((r^), 
respectively. The corre- 
^ 0.595 is obtained from 
as marked by the circle. 



appear to be rather unnatural. It is not possible in the 
ghost-antighost symmetric formulation. 

Another important difference between 1^1 and 1% for 
the tree-level vertex, as seen in Fig. [|and Eq. (|8^), is the 

observation that 1^1 (k) has a pole at k = 1/2. The gluon 
propagator then necessarily vanishes in the infrared: If 
it was to approach a constant, one had to have k = 1/2. 
In this case, however, its constant limit was proportional 
to 1/7^2 (k) which vanishes for k — ► 1/2. One thus ob- 
tains the strict lower bound 1/2 < k, for the tree-level 
case. Zwanziger's horizon condition is then satisfied. The 
apparently infrared finite extrapolations from lattice cal- 
culations are an open question still, however. 



B. Infrared transversality 

So much for the tree-level vertex. Before we discuss 
other, more gen eral p ossibilities, in particular, the cases 
(i) - (iii) in Sec. Ill B| , we give a more detailed discussion 
of the transversality issue in this subsection. 

For the "^-contracted infrared contribution of the 
ghost-loop in the gluon DSE, Eq. (|3l|), we must specify a 
form for the second, the longitudinal -B-structure of the 
vertex, in addition. The fact that the leading infrared 
behavior of the propagators should not depend on the 
choice of studying Zp or Zr can be used to construct an 



infrared form of B(k ,p , q ) analogous A %r in (p6|). 

First, we express the integrand in (|3l]) for Zp in terms 
of that for Zp in (30) plus a correction term which, for 
a given Ansatz A lr , is required to vanish (at least upon 
integration). For the four terms in the curly brackets in 



(pl|), this leads to 

{•••} = qV(k)qA(k 2 ;q 2 ,p 2 ) - (D - 1) X (83) 
(f)k 2 + qk) ^(l + g^fcW,^) - B(fcW,p 2 )) • 

The first term on the r.h.s herein reproduces the result 
for I^J (k, I, n) from (|74|). With tree-level settings A = 1, 
B = f\ one readily verifies Eq. (82). 

The difference between the infrared integrals, A( D > = 



t(D) 



— I^PK can then be written, 



2 D w D / 2 
4 l 



p_ 

k 2 



(2tt 

2 „2 



D 



1+ft 



(1 



El 

k 2 



Si' 
k 2 



A' tr - 2B ir . (84) 



The order of the arguments in A %r and B %T is the same 
here as that in ( |S3| ) above. The ghost-loop integra- 
tion projects on terms that are overall symmetric in 
p 2 <-> q 2 . The antisymmetric ones vanish upon inte- 
gration. We use the symmetry (SI) of A and the de- 
composition B = B + + B^. into (anti)symmetric parts 
B±(x; y, z) — ±B±(x; z, y) for B again, and sort out the 
symmetric part of the integrand in (p4[). It vanishes, if 

(fj- V )A^(k 2 ]q 2 ,p 2 ) = (f 1 -r ] )2B^(k 2 ;q 2 1 p 2 ) , and 
p 2 -q 2 



2k 2 



A-(k 2 ;q 2 ,p 2 ) = BZ{k 2 -q 2 ^). 



(85) 



The first condition shows that 2B^ = A %r , for i) ^ r\. 
In the symmetric case, 77 = 77 = 1/2, no such restric- 
tion is implied her e, bu t the same relation is then given 
by (S2), see Sec. |lA| above. Therefore, Eqs. (§2|) and 
( p3| ) of Sec. II A follow as sufficient conditions for infrared 
transversality independent of rj. Without turning them 
into necessary cond ition s this is a rather trivial result. As 
we showed in Sec. [I A , the conditions in ( ^2|) and (|23| ) 
imply that the vertex itself is transverse in the infrared. 
This is always a possibility to warrant infrared transver- 
sality, however. The necessary condition is = 0. 

The point here is to demonstrate that, apart from pos- 
sible accidental cancellations, e.g., with tuning 77 = 1.16 
for the tree-level case of the last subsection, for general 
77, there are really no possibilities left other than in- 
frared transversality of the vertex itself. In particular, we 
should be allowed to choose 77 such as rj = 1 for the stan- 
dard linear-covariant or 77 = 1/2 for the ghost-antighost 
symmetric case. We might then further say that we are 
not interested in contributions to the vertex which them- 



selves vanish upon integration in Eq. (84) for A( D \ Up 
to such irrelevant contributions, which neither contribute 
to the gluon nor the ghost DSE, ( p2| ) and ( p3| ) and thus 
the transversality of the vertex in the infrared are also 
necessary conditions for the infrared transversality of the 
gluon correlations in Landau gauge independent of 77. 
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FIG. 4: Roots without symmetrization (left), real roots with additional symmetrization (right). 



C. Bounds on the infrared exponent 

We now go back to the general results given in Eqs. ( |7q ) 
and (|7r|). With these results for the necessary in- 
frared integrals it is little effort to explore infrared 
forms A lr (k 2 ;p 2 , q 2 ) other than A lr — const, for the A- 
structure in the ghost-gluon vertex. First, in 4 dimen- 
sions, the selfconsistency condition in Eq. ( |72| ) for these 
integrals leads to, 

(l+n+K-l)(l+n+K)(l+n+K+l)(l+n+K+2) = 

-3(n + 2/t)(n+2K-l)(n+2K-2)(n + 2K-3) . (86) 

This then corresponds to the Ansatz for the vertex as 
given in (|56|), (57) with only condition (Nl) being im- 
plemented at this stage. It is the starting point for the 
discussion of the t hree special cases introduced in Eqs. 
(H) " © of Sec. |IIIB| : 

Case (i). Here, we simply need to set n+l = 0. We then 
obtain form ( p6| ) for the possible solutions to Eq. (|7£ 
which here reads, 

I^ ] (K,-n,n) = I ( zl(n,-n,n) , (87) 

the quartic equation to, e.g., determine n(n), 

(k-1)k(k+1)(k+2) = -3(u+2k) ■ ■ ■ (n+2 K -3) . (88) 

The four real roots to this equation are shown for n in 
[0, 1] in Fig. | on the left. 

The top branch, with 1 < n < 3, leads to negative 
Iq and 1$. Both propagators, Z{k 2 ) and G(fc 2 ), then 
necessarily had zeros at some finite fc 2 , and it resulted 
that a c < (equality at the bounds n = 3, 1 for K = 0, 1). 
We therefore rule out this solution as being unphysical. 

Also, we are particularly interested in solutions with at 
most weak singularities in vertex functions and require 



therefore \n\ < 1. This is the case for the next branch 
with < n < 2 provided 0.4222 < k. This branch will no 
longer exist after symmetrization w.r.t. ghost- antighost 
momenta (see Case (ii) below), however. 

For the bottom branch with < n < —2, the re- 
striction to weakly singular vertices, — 1 < n leads to 
k < 0.4767. With this solution, it is therefore impossible 
to obtain the mass gap in transverse gluon correlations, 
for which 0.5 < k. Furthermore, it does not survive the 
symmetrization in (ii) either. 

For the branch that includes the tree-level result, that 
with n(K\) — 0, the critical coupling from (f70|), with 
D = 4 and N c = 3, is shown as a function of k in Fig. ^ 
(dashed line). Its maximum a™ ax 2.9798 occurs at 
k w 0.6174. It is thus slightly larger than the tree-level 
value a c (Ki) » 2.9717 given in Eq. (|l|). 

Case (ii). This case corresponds to a sum of two terms, 
one with m = and I = —n as in Case (i) above, and 
the other with / = and m — —n. Since the ghost-loop 
contribution Iz P is manifestly symmetric inp 2 <-> q 2 , this 
symmetrization only affects the infrared contribution Iq 
to the ghost DSE, and we can therefore write, 

1 



I G {n,n) 
I Zp (t, n) 



4 4) ( K ,-n,n) + 7^,0,71) 



The selfconsistency condition Iq(k,ti) = Iz P {n,n) 
now be used to obtain, 

(k-1) ■ ■ ■ (k+2) (n+K-1) ■ • • (n+K + 2) \ 



(89) 

(90) 
can 



(k-1) • • • (k+2) + (n+K-1) ■ ■ ■ (n + K+2) 
--(n + 2n)(n+2n-l)(n+2n-2)(n+2n- 



3). (91) 



This equation has eight roots n(n) which, in general, 
come in complex conjugate pairs. The real roots for k 
in [0, 1] are shown in Fig. |] on the right. 
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FIG. 5: The critical value of the running coupling over the 
infrared exponent re for Case (i) (dashed) and Case (ii) (solid). 
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FIG. 6: Solution n(re) for Case (m) (dashed) compared to the 
tree- level branch of Case (ii) (solid) in the range 1/2 < re < 1, 
both with n(rei) = at rei = (93 - \/T20T)/98. 



One discovers that two out of the originally four 
branches, from Case (i) above, are almost unaffected by 
the symmetrization employed here: These are the un- 
physical one at the top, with 1 < n < 3 and much the 
same a c < 0, and the one connected to the tree-level re- 
sult, with ti(k\ ) = for ki = (93-Vl201)/98. Superim- 
posing both results for the latter, with and without the 
symmetrization, the corresponding solutions n(n) turn 
out to be almost indistinguishable in the whole range 
< At < 1 on the scales of our plots. A little bit more 
appreciable, but not very significant still, are the differ- 
ences in the corresponding values for a c as compared to 
each other in Fig. |5|. For the symmetric solution (solid 
line) the maximum occurs at the value Ati for the tree- 
level vertex, a c nax = a c («i), at whic h both solutions in- 
tersect. As explained in Sec. Ill B , if we furthermore 
require the vertex-structure A to have no infrared di- 
vergences associated with the ghost momenta, we must 
have n < for the ghost-antighost symmetric vertex, in 
addition. Therefore, we can find physically acceptable 
solutions in the range 



«i < At < 1 , and < a c (/t) < a c (Ati) 



(92) 



correspondingly, with the values Ati ~ 0.59535 and 
Qe(Ati) ~ 2.9717 of the bounds as obtained from Eqs. j80| ) 
and ( j8l| ) for the tree-level vertex, respectively. 

For completeness we mention that the new branch for 
the symmterized vertex (ii) with —1.2 < n < 0, with 
possibly interesting solutions — 1 < n for At < ti> leads 
to a c < throughout. And this is also the case for the 
bottom branch with n < — 2. 

In the range of particular interest, 1/2 < n < 1, we are 
thus left with the branch of solutions including u(ki) = 
as the only one with a c > after the symmetrization in 
(ii). At the same time, it seems quite encouraging that 
this branch, the only physically relevant one, is practi- 
cally unaffected by the symmetrization. 

Case (Hi). In this example, the infrared behavior of 
the ghost-gluon vertex as given in Eq. (|6l|) is such 



that it again satisfies the conditions (Nl), from non- 
renormalization, and (SI), from ghost-antighost symme- 
try. In the limit where one ghost momentum vanishes, 
one now has A(q 2 ;q 2 ,0) = 0, however. We will see that 
this has no dramatic consequences either, on the physi- 
cally interesting solutions found in the range (92) above. 

Here, by the same arguments as in Case (n), we can 
now express the leading infrared integrals in both DSEs, 

I G ( K ,n) = 4 4) («,-n,n) + 4 4) (AC,0,n)-4 4) (A t ,0,0) , 
I Zp (At, n) = 2lf p (/t, -n, n) - lf p (At, 0, 0) . (93) 

Due to the (I = n = 0) contributions herein, which arise 
from the tree-level term (with negative sign) in (|6l|), it 
is generally no longer possible to derive the solutions to 
Ig{k, n) = Iz p (k, n) as the roots of simple polynomials. 
Searching the physically interesting range of parameters 
numerically, starting from the known solution for n = 0, 
k — Ati, we obtain the dashed curve for n(At) as compared 
to the corresponding branch for Case (ii) in Fig. ^. Again 
requiring n < to avoid infrared divergent ghost legs, we 
find that the solutions in the two cases are remarkably 
close to each other with — k < n < for Ati < At < 1. 
And again we find that n — > — k in the limit At — » 1 in 
which a c — > in all three cases, however. 

One might think that another solution with n = 
exists for At = 1/2 and the Case (Hi) vertex. Since that 
could have important implications, we note here that this 
is actually not the case. We know that Iz P in (^3|) re- 
duces to the form for the tree-level vertex at n — 0, which 
does not lead to a solution in 4 dimensions, c./., Fig. ||. 
The Fact that the dashed line in Fig. || appears to ap- 
proach n — > + for At — > (l/2) + is explained as follows: 
For sufficiently small n = e, we find from (193) that 



(At,e) 



1 



1 



1 



At-(l-e)/2 2At-l/2 



for At 



1 

2 ' 



is dominated by two nearby Poles with opposite signs. 
Therefore, for arbitrarily small but finite e > the Pole 
at At = 1/2 in the contribution from the last term, the 
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FIG. 7: Iz P (K,n) (dashed) for several several small values 
n = {0.01, 0.006, 0.003} and n = over the infrared exponent 
k in Case (Hi). The intersection with Ig(k, 0) (solid) near 
K = 1/2 disappears for n = 0. 
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FIG. 8: The value of a c in the range 1/2 < k < 1 for all 
three Cases, (i) (long dashed), (ii) (solid), and (Hi) (short 
dashed) . 



negative of that in the tree-level vertex case, will always 
lead to an intersection, just above n = 1/2, with Ig(k, 0) 
which approaches a constant (corresponding to the value 
47r/(37 G (0.5, 0)) w 2.62 in Fig. |). For n = e = on the 
other hand, both Poles coincide and their residues sum 
up to that of the tree-level case which is now positive. 
And thus, the intersection point then disappears. This is 
confirmed also numerically and demonstrated in Fig. 0. 

Of course, if we relax the condition n < 0, we can 
have selfconsistent solutions also for infrared exponents 
k < «i, including those for k = 1/2 in Cases (i) and 
(ii). At the same time, this leads to a singularity in 



A(k 2 ;p 2 ,q 2 ) as q 



0. 



Negative n on the other hand lead to an infrared di- 
vergence associated with the gluon leg. However, as long 
as — n < k this is over-compensated by the gluon propa- 
gator attached to that leg, because Z(k 2 ) ~ (k 2 ) 2K . For 
n = —K an effective massless particle pole would be left in 
a gluon exchange between two vertices, G ll D fliV (k)G l , ~ 
1/fc 2 . In all solutions we report here, n — > —k for 
ct c ~~ * (both from above). Therefore, this limit can- 
not be reached, since then, at least one of the leading 
infrared coefficients do or eo in the propagators vanishes 
which contradicts the assumptions, c./., Eqs. (|53|) (|55|). 

No such compensation occurs for divergences associ- 
ated with ghost legs. The ghost correlations are them- 
selves infrared enhanced. This infrared enhancement will 
persist for the ghost correlations between their vertices, 
if < k — m — I < 1. Above the upper bound, the in- 
frared divergences become too severe for the description 
in terms of local fields. For the cases with the ghost- 
antighost symmetry (SI) of Landau gauge, we obtain 
from this restriction the upper bound n < (1 — k)/2 < 
0.25 for k < 1/2. This, however, leaves just enough room 
that it alone does not rule out the solutions with positive 
n found for 1/2 < k < K\ in case (ii), as seen in Fig. ^| 



V. SUMMARY AND CONCLUSIONS 

The Dyson-Schwinger equations of standard Faddeev- 
Popov theory in Landau gauge, when supplemented by 
additional boundary conditions, can be derived as an 
approximation to the time-independent diffusion equa- 
tion of stochastic quantization which is valid non- 
perturbatively |j35f . The non-conservative part of the 
drift force that is neglected in this approximation can- 
not be described by local interactions. The effect of this 
part will have to be investigated in future. It may well 
be responsible for the kind of "Gribov noise" observed in 
lattice calculations. 

Here, we studied a slightly more general definition of 
the Landau gauge as limit of a wider class including 
non-linear covariant gauges |52[] , This limit is controlled 
by an additional free parameter 77 in the tree-level ver- 
tex (with X] = 1 in Faddeev-Popov theory). In particu- 
lar, we find that non-renormalization of the vertex in a 
symmetric subtraction scheme and infrared transversal- 
ity of the gluon propagator in Landau gauge can only go 
together with the manifestly ghost-antighost symmetric 
choice 77 = 1/2. In the light of the recent progress con- 
necting the linear-covariant gauge with time- independent 
stochastic quantization, the ghost-antighost symmetric 
Curci-Ferrari gauges might therefore also deserve to be 
reconsidered for similar connections. 

Optimistically assuming that perfect sense can be 
made of Dyson-Schwinger equations non-perturbatively 
some day, we studied the infrared critical exponent and 
coupling for gluon and ghost propagation in Landau 
gauge in quite some generality. We gave two reasons 
for assuming ghost dominance, the Kugo-Ojima criterion 
for confinement and the horizon condition to restrict the 
measure to the first Gribov region, and implemented this 
as a boundary condition in our infrared asymptotic dis- 
cussion of the DSE solutions. Central to an understand- 
ing of the infrared exponents for gluon and ghost propa- 
gation in Landau gauge then is a knowlegde of their ver- 
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tex. Assuming it has a regular infrared limit, we obtain 
k ~ 0.595. For the ghost-antighost symmetric vertices, 
this value maximizes the critical coupling a c (n), yielding 
a max w 2.97 ; as summarized once more for k between 0.5 G f j^ e f 54 



APPENDIX A: Ghost RG equation revisited 



and 1 in Fig. pi For larger k the vertex acquires an in- 
frared singularity in the gluon momentum, smaller ones 
imply infrared singular ghost legs. 

Quite encouragingly, numerical solutions to truncated 
Dyson-Schwinger equations have recently been obtained 
with the infrared behavior of the regular vertex, as 
obtained here, for the whole momentum range up to 
the perturbative ultraviolet regime and without one- 
dimensional approximation in Ref. [ p5| . 

An important detail in all our considerations is the 
non-renormalization of the ghost-gluon vertex in Landau 
gauge. Derived from standard Slavnov- Taylor identities 
it is one of the arguments that hold at all orders of per- 
turbation theory. That it is also true non-perturbatively, 
however, is another additional assumption. It is there- 
fore quite important and interesting that this has been 
assessed and verified within the numerical errors in cal- 
culations using Landau gauge on the lattice p6| . Calcu- 
lating both propagators simultaneously, this study fur- 
thermore appears to confirm a unique exponent for the 
combined infrared behavior of gluons and ghosts for the 
first time in a lattice calculation implementing the Lan- 
dau gauge condition. Also, for SU (2) this study reports 
preliminary values of a c that are fully consistent with 
the results obtained here somewhere near the maximum 
a rnax _ (4^/^)0.71 ~ 4.5 for N c = 2. It will be very 
interesting to see the final errors, so that we will be in the 
fortunate position to restrict further the range of both, a c 
and k. At the moment, the combined evidence seems to 
indicate that the result will be somewhere in the range 
around k — 0.5 and the maximum near k — 0.6. Un- 
fortunately, with this conclusion the question about an 
infrared vanishing versus finite gluon propagator must 
therefore remain open, for the time being. 
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consider formula (Bl) as new in any sense. This is not 
at all the case. The 2-line derivation in (B2), (B3) below 
is given for the convenience to the reader. 

He furthermore points out that relation ( B13| ) follows 
from Eq. (11) listed on page 534 in the tables of Ref. Jro"| . 
We gratefully acknowledge this information. 



We repeat the renormalization group (RG) analysis 



for the ghost propagator with some minor 
corrections. These corrections do not affect the main 
conclusions of Ref. [Q as far as we judge. The correct 
versions, in particular, of Eqs. (8) to (11) in Ref. Q, 
are necessary, however, in order to establish the equiv- 
alence of their asymptotic infrared expansion and the 
one adopted via Eq. ( (46| ) herein, which is a minor vari- 
ation of the expansion techniques developed previously 

@ f§ g II g |g, |§. 

First, recall the RG equation for the ghost propagator 
G(k 2 ) = G(k 2 , [i 2 ), in this form also given in (6) of Q, 



2 lG (g) \G(k\^)=0 



(Al) 



The formal solution to this equation is given by Eq. (|41 
With the infrared Ansatz in the form of Eq. (5) in |54 



N 



G{k 2 ^ 2 )^Y, d '^){- 2 



(A2) 



here denoting the coefficients of [54] as the primed ones, 

(4£), we first obtain, 



d' n , to distinguish from those in 



/?(<?) 



9< 
dg 



d' n Hk 2 /^ 



d5 n 
dg 



2d' n {8 n 



•7c(ff)) =0, 
(A3) 



at variance with Eq. (8) of Ref. |54|] in two minor ways (by 
the factor of 2 and the sign of the 7c(g)-term). Never- 
theless, with their conclusion that therefore dS n /dg = 0, 
we find for the coefficients, 



dd' n 2(5 n + lG (g)) 
dg P(g) 

the general solution of which takes the form, 



(A4) 



[ a 2(6 n + lG (l)) M .... 
d n (g) cx exp J — '-dl . (A5) 

This, however, is incompatible with Eq. (9) of Ref. |p3, 



2(S n +y G ) 

d n (g) = const, g 2 ~<g+~ia . 



(A6) 



In particular, since we just noted that the exponents 
5 n are (/-independent, they are eithe r ze ro or o ne w ould 
need j3{g) = —const, x g to obtain (|A6[ ) from (|A5j ). By 

virtue of Eq. @, 0(g) = -g{2 lG + lA ) g, this 

would imply that 2"f G + jA is (/-independent. If we would 
then conclude in addition that both, j G and 7^, a re g - 
independent, only then we would obtain ( A6) from ( |Aq ). 

Note that such a behavior, (3(g) ~ —g, though in prin- 
ciple possible in the infrared, would not lead to a fixed 
point and thus contradict the other results of p3| as we 
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discussed in Sec. Ill A above. Not restricted to such a 



specific behavior, here we go back to the general form 
of the d' n in Eq. (A5). First, remember that g = go for 



/i = a . In this case, the exponential factor in our expan- 



sion ( fig ) becomes unity, and ( J4q ) and (A2) a gree . Thus, 
d n = d' n (go), and we can split the solution to (A4) for d' n 
with this initial condition into factors as follows, 



d 'n(g) = d n exp{25 n 



£|y }exp{2 £ 



9 i^idD 



(A7) 



The first exponential factor herein, with ([44]) , is equal to 
(^ 2 /o-) 5 ", and can be used to replace — > a in flA^). 
The last exponential in (A7) is the same overall factor 
determining the (g, /^-depe nden ce of the ghost propaga- 
tor as in (|46|). Substituting (A7) into the expansion ( |A2| ) 
of Ref. [ML one obtains, 



G(k 2 ,fi 2 ) ~ exp{2 



9 1g{1) ... , fk^ 5 



dl } d n 



(A8) 



0(0 

which agrees with the RG invariant expansion of (E(t) 



APPENDIX B: Two ways to do the D-dimensional 
integrals 

The basic formula we employ for the infrared analysis 
in Sec. Ill involves D-dimensional integrals of the follow- 
ing form which converge for Re(a) > 0, Re(/3 — a) > 0, 
Re(/3) < D/2, 



(2ir) D \q 



D/2 



(Bl) 



1 T(a)T{D/2 - f3)T{P - a) 

2 D w D / 2 T(fi)Y{Dj2 - a)T{D/2 +a-0) ' 

where p = k ± q and an explicit factor (k 2 ) ,3 ~ a was in- 
troduced to render the integral dimensionless. This is a 
textbook formula, c./., Eq. (2.5.178) in 163). For a simple 



derivation with t he g eneral exponents |64| one observes 
that the l.h.s. of (Bl) is a convolution integral which re- 
duces to an ordinary product upon Fourier transforma- 
tion. Using 



1 r(D/2- 7 ) 



(g2)7 477T D / 2 r( 7 ) 



dPx {x 2 y~ D/2 e~ iqx , (B2) 



for the two factors in the convolution (with the power 7 
given by D/2 — a and /3, respectively), one thus obtains, 



d D q fl\ D/2 
(2ir) D \q 2 ) \k 2 



4 a -P T(a)T(D/2- f3) 
(2tt) d T(p)T(D/2-a) J (x 2 ) D / 2 



(k 2 x 



2 2\f3-a 



(B3) 



-ikx 



Eq. (Bl) then follows from a further application of the 
Fourier transform (B2) herein, now with 7 = (3 — a. 

Alternatively, we can do the integral in (Bl) which 
we denote henceforth by fu{ot,f3) in a straightforward 
though less elegant way by first performing all but one 
of the angles of the polar coordinates in _D-dimensional 
momentum space, 



fD(a,0) 



K(D) D f dq 



(2ir) D 2 J q 2 \ k 2 



(B4) 



1 



B((U-l)/2,l/2) 



dz{\- z z ) 



2\(£>-3)/2 / ^_ 



2\ P 



P 



where K(D) = 2tx d I 2 / \DY{D /2)) is the volume of the 

I 2 

Euler Beta function is given by 



D-dimensional unit ball, p 2 = q 2 + k 2 — 2kqz, and the 



£((£>- 1)/2, 1/2) = f rfti- 1 / 2 (l - £ )(^-3)/2 
Jo 



dz(l 



2x(D-3)/2 



(B5) 



For the azimuthal integration in (B4), the formula 2. in 
3.665 of Ref. |65| can be used to obtain, 



f D (a,/3) 



K(D) D 
(2^2 



1 



2 D n D / 2 T(D/2) 



(B6) 



2 D ir D / 2 T(D/2) \a 



i 3 F 2 (/3,/3™ + l,a;~,a+l;l) + -g^ 3 F 2 ^,P-j + l ) /3-a;^,p-a+l;lj S j , (B7) 



r 



where for the last step the following integration formula final result in (B7) 
for the generalized hypergeometric functions led to the 



(B8) 



/ dxx 1 pF q (a;h;x) — 
Jo 

— j— p+1 F 9+1 ({a, 7 +l};{b, 7 +2};l), 
7+1 
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which is most easily derived for 7 > — 1 (see, e.g., the 
appendices of Refs. jh| |47| ) from the power series ex- 
pansion for the generalized hypergeometric functions, 



p F q (a;b;z) = ^ 



n=0 



by noting the relation, 
1 



(fll)n ' ' ' { a p)n 

(bi) n ■ ■ ■ {bq) n n\ 



1 (7 + l)r 



7 + n + 1 7 + 1 (7 + 2) 7 
for the Pochhammcr symbols, 

T(a + n) 



{a) r 



T(a) 



(B9) 



(BIO) 



(Bll) 



For other general properties and a variety of relations 
amongst the different hypergeometric functions, e.g., see 
Refs. §,§§§. A well-known one for example 
expresses the Gauss series as a ratio of Gamma functions, 



2 Fx(a,b; c; 1) 



T(c)T(c- 



b) 



T(c-a)T(c-b) 



(B12) 



Many more relations of this kind, including less known 
ones, are listed in the tables on hypergeometric functions 
of Ref. [0 . We originally thought it might be interesting 
to note that the two ways to calculate fo{a,P) leading 
to the r.h.s. of Eqs. (Bl) and (B7), respectively, allow to 
devise some of these additional relations. For example, 
by simple renamings, a comparison of the r.h.s. in (Bl) 
and (B7) yields, 



- 3 F 2 (a,a+b,a+b-c+l;c,a+l;l) + \ 3 F 2 (b, a+b, a+b-c+1; c, 6+1; 1) 



T(a)T(b) T(c)T(c-a-b) 
T(a + b) T{c-a)T(c-b) 
B{a,b) 2 F 1 (a,b;c;l) . (B13) 



This formula follows with replacing b — > a + b, 1 — c ^ 
c — a — b, and d — > c upon rearrangement from Eq. (11) on 
p. 534 in |7Q], and our presentation here seems obsolete 
now 72 . At least, the equivalence of the results from 
the infrared a naly sis of DSEs in Ref. fig] , to the expres- 
sions in Eqs. ([77]) a nd (|82] ) with -q — for the tree-level 
vertex case of Sec. IV A is explicitly established in this 



way. The second procedure to calculate integrals such as 
/i)(a,/3) was thereby used in Ref. Each of the re- 

sults therein are readily expressed in terms of one simple 
ratio of Gamma functions when using the relations pre- 
sented in this appendix. Though equivalent to Eq. (|B7]) 
of course, use of Eq. (Bl) thereby is far more convenient 
for all practical purposes. 
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